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Rickettsia  prowazekii,  the  causative  agent  of  epidemic  typhus,  has  been  responsible  for  millions  of  human 
deaths.  Madrid  E  is  an  attenuated  strain  of  R.  prowazekii,  while  Breinl  is  a  virulent  strain.  The  genomic  DNA 
sequence  of  Madrid  E  has  recently  been  published.  To  study  the  genomic  variations  between  Madrid  E 
(reference)  and  Breinl  (test)  DNAs,  cohybridization  experiments  were  performed  on  a  DNA  microarray 
containing  all  834  protein-coding  genes  of  Madrid  E.  Of  the  834  genes  assessed,  24  genes  showed  1.5-  to  2.0-fold 
increases  in  hybridization  signals  in  Breinl  DNA  compared  to  Madrid  E  DNA,  indicating  the  presence  of 
genomic  variations  in  ~3%  of  the  total  genes.  Eighteen  of  these  24  genes  are  predicted  to  be  involved  in 
different  functions.  Southern  blot  analysis  of  five  genes,  virB4,ftsK,  rfbE,  IpxA,  and  rpoH,  suggested  the  presence 
of  an  additional  paralog(s)  in  Breinl,  which  might  be  related  to  the  observed  increase  in  hybridization  signals. 

Studies  by  real-time  reverse  transcription-PCR  revealed  an  increase  in  expression  of  the  above-mentioned  five 
genes  and  five  other  genes.  In  addition  to  the  elevated  hybridization  signals  of  24  genes  observed  in  the  Breinl 
strain,  one  gene  ( rp084 )  showed  only  1/10  the  hybridization  signal  of  Madrid  E.  Further  analysis  of  this  gene 
by  PCR  and  sequencing  revealed  a  large  deletion  flanking  the  whole  rp084  gene  and  part  of  the  rp083  gene  in 
the  virulent  Breinl  strain.  The  results  of  this  first  rickettsial  DNA  microarray  may  provide  some  important 
information  for  the  elucidation  of  pathogenic  mechanisms  of  R.  prowazekii. 


Rickettsiae  are  obligate  intracellular  gram-negative  bacteria 
that  belong  to  the  alpha  subdivision  of  Proteobacteria  (48). 
Rickettsial  diseases  are  widely  distributed  throughout  the 
world  and  are  characterized  by  sudden  onset  of  acute  fever, 
skin  eschars,  headache,  and,  in  cases  of  the  fulminant  disease, 
disseminated  intravascular  coagulation  that  leads  to  failure  of 
the  cardiovascular  system.  Diagnosis  of  rickettsial  diseases  is 
difficult  due  to  the  lack  of  specific  signs  and  symptoms  and 
the  inability  to  grow  the  bacterium  on  axenic  microbiological 
media  for  confirmation  (32).  There  are  no  licensed  vaccines 
available  for  the  prevention  of  rickettsial  diseases  in  the 
United  States.  Moreover,  three  rickettsial  agents,  Rickettsia 
prowazekii,  Rickettsia  rickettsii,  and  Coxiella  burnetii,  have  been 
listed  as  biothreat  agents  by  the  Centers  for  Disease  Control 
and  Prevention  in  the  United  States. 

R.  prowazekii  is  the  etiologic  agent  of  epidemic  typhus  in 
humans.  Its  presence  is  commonly  associated  with  war,  famine, 
and  social  disruption.  Typhus  killed  >3  million  people  and 
infected  in  excess  of  20  million  people  during  World  War  I, 
and  it  killed  several  million  more  in  World  War  II  (20).  As 
recently  as  1997,  there  was  an  epidemic  of  typhus  in  Burundi 
that  infected  well  over  100,000  people  (31).  Other  outbreaks  in 
Eastern  Europe,  Russia,  Ethiopia,  and  South  America  further 


*  Corresponding  author.  Mailing  address:  Rickettsial  Diseases  De¬ 
partment,  Infectious  Diseases  Directorate,  Naval  Medical  Research 
Center,  Silver  Spring,  MD  20910.  Phone:  (301)  319-7438.  Fax:  (301) 
319-7460.  E-mail:  chingw@nmrc.navy.mil. 


support  the  view  that  epidemic  typhus  is  a  reemerging  infec¬ 
tious  disease  worldwide  (20,  21).  Present-day  threats  posed  by 
typhus  bacteria  include  natural  exposure  to  antibiotic-resistant 
strains  and  the  use  of  laboratory-manipulated  strains  as  bio¬ 
logical  warfare  agents  (28).  The  World  Health  Organization 
has  estimated  that  104,000  casualties  (19,000  dead  and  85,000 
incapacitated)  would  arise  following  a  hypothetical  air  deploy¬ 
ment  of  50  kg  of  typhus  agent  (50). 

The  pathogenic  mechanism  of  rickettsiae  is  unclear  but  may 
involve  potential  virulence  factors,  such  as  phospholipases  (37, 
46,  49),  antigenic  surface  proteins  (10,  19,  24,  46),  hemolysins 
(29),  type  IV  secretion  systems  (2,  25,  26),  and  the  recently 
described  nudix  hydrolase  (18).  Significant  biological  differ¬ 
ences  exist  between  the  virulent  Breinl  strain  of  R.  prowazekii 
and  the  attenuated  avirulent  Madrid  E  strain  (42).  The  Breinl 
strain  grows  well  in  macrophage -like  cells,  while  the  Madrid  E 
strain  grows  poorly.  The  Breinl  strain  has  a  greater  cytotoxic 
effect  on  macrophage -like  cells  treated  with  alpha/beta  inter¬ 
feron  (IFN-a/(3)  or  IFN-y  than  the  Madrid  E  strain,  and  the 
Breinl  strain  is  less  sensitive  to  IFN-y  treatment  than  the 
Madrid  E  strain  in  L929  cells  (42).  However,  the  genes  respon¬ 
sible  for  the  variation  in  virulence  between  the  two  strains  of  R. 
prowazekii  have  not  been  well  characterized.  The  published 
genomic  sequence  of  the  Madrid  E  strain  (2)  provides  a  foun¬ 
dation  for  comparative  genomic  studies  of  R.  prowazekii.  The 
circular  genome  consists  of  1.1  Mbp  containing  834  protein¬ 
coding  genes,  12  pseudogenes,  and  38  RNA  genes.  About  24% 
of  the  genome  belongs  to  the  noncoding  region,  which  repre- 
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sents  the  highest  proportion  of  noncoding  DNA  detected  in 
any  microbial  genome  sequenced  so  far  (2).  The  comparisons 
of  functional  profiles  of  Madrid  E  genes  with  those  of  mito¬ 
chondrial  genes  revealed  many  similarities  in  ATP  production 
and  transport  systems  (2). 

The  purpose  of  these  studies  was  to  compare  the  genomic 
content  of  the  avirulent  Madrid  E  strain  with  that  of  the  vir¬ 
ulent  Breinl  strain.  We  report  the  construction  of  an  R. 
prowazekii  DNA  microarray  with  834  PCR-amplified  gene 
products  from  Madrid  E  and  the  comparative  genomic  studies 
with  DNAs  from  Madrid  E  and  Breinl  using  this  microarray. 
The  microarray  results  led  to  the  discovery  of  a  large  deletion 
including  the  rp084  gene  in  the  Breinl  strain.  Southern  blot 
analysis  of  five  dissimilar  genes  revealed  the  possible  presence 
of  additional  paralogs,  and  the  increased  expression  of  10 
virulence-related  genes  was  confirmed  by  real-time  reverse 
transcription  (RT)-PCR.  Elucidation  of  the  roles  of  these 
genes  could  lead  to  a  better  understanding  of  the  mechanisms 
for  the  pathogenesis  of  R.  prowazekii.  In  addition,  insights 
gained  from  these  studies  may  have  implications  for  the  devel¬ 
opment  of  new  diagnostic  strategies,  treatment  targets,  and 
vaccine  candidates  for  epidemic  typhus. 

MATERIALS  AND  METHODS 

Rickettsial  strains  and  purifications.  R.  prowazekii  strain  Breinl  was  originally 
isolated  from  a  patient  in  Warsaw,  Poland,  in  1922  and  was  brought  to  the 
United  States  in  1928.  This  strain  was  passaged  multiple  times  and  used  as  a 
representative  strain  of  epidemic  typhus.  The  strain  used  in  this  study  was  in  its 
~  157th  passage.  Madrid  E  was  the  egg  passage-attenuated  strain  of  the  Madrid 
strain,  which  was  originally  isolated  from  a  typhus  patient  in  Madrid  (2).  The 
Madrid  E  strain  used  here  was  in  its  —281st  passage.  The  exact  passage  histories 
of  the  two  strains  are  unknown  because  of  the  lack  of  data  collected  early  in  the 
cultivation  of  the  strains. 

Six-day  pathogen-free  chicken  embryos  were  inoculated  with  yolk  sac  suspen¬ 
sions  from  ampoules  and  incubated  at  35°C.  The  infected  yolk  sacs  were  pooled, 
shaken  with  glass  beads  and  further  macerated  by  pipetting,  diluted  with  brain 
heart  infusion,  and  centrifuged  at  10,000  X  g  for  30  min.  The  pellets  were  washed 
twice  with  buffer  K36  (0.1  M  KC1,  0.015  M  NaCl,  0.05  M  potassium  phosphate, 
pH  7.0)  by  resuspension  and  sedimentation  at  17,300  X  g.  The  pellet  was  sus¬ 
pended  in  one-half  volume  of  the  original  brain  heart  infusion,  and  5  ml  of  this 
suspension  was  layered  over  20  ml  of  10%  bovine  albumin  in  buffer  K36.  After 
centrifugation  at  480  X  g  for  20  min,  small  sediments  were  discarded  and  the 
supernatants  were  diluted  with  3  volumes  of  K36  and  centrifuged  at  10,000  X  g 
for  30  min.  The  washed  pellets  were  resuspended  in  K36  and  subjected  to  two 
consecutive  cycles  of  Renografin  density  gradient  purification  (13).  In  both  cases, 
2.5  ml  of  rickettsial  suspension  was  layered  over  the  following  gradients,  (i)  The 
first  was  a  30  to  45%  linear  gradient  of  Renografin  for  centrifugation  at  25,000 
rpm  (SW27  rotor)  for  1  h.  After  centrifugation,  the  areas  of  the  gradient  above 
and  below  the  rickettsial  band  were  removed  with  a  Pasteur  pipette.  Then,  the 
rickettsial  band  was  drawn  into  a  syringe  through  a  14-gauge  cannula,  diluted  10 
times  with  K36,  and  centrifuged  at  17,300  X  g  for  10  min.  The  pellets  were 
suspended  in  K36.  (ii)  The  second  was  a  20  to  45%  linear  gradient  of  Renografin 
with  centrifugation  at  25,000  rpm  (SW25.1  rotor).  The  suspension  from  the  first 
gradient  centrifugation  was  processed  through  the  second  gradient  centrifuga¬ 
tion.  The  pure  rickettsial  band  was  drawn,  and  the  final  pellet  was  suspended  in 
K36.  Live  rickettsiae  purified  by  this  method  were  essentially  free  of  host  mate¬ 
rial. 

Construction  of  microarray  slides,  (i)  DNA  extraction  and  amplification. 

Renografin-purified  R.  prowazekii  Madrid  E  was  resuspended  in  ATL  lysis  buffer 
(Qiagen,  Valencia,  Calif.).  After  proteinase  K  digestion  and  RNase  A  treatment 
according  to  the  instructions  of  the  manufacturer,  DNA  was  extracted  three 
times  with  phenol-chloroform-isoamyl  alcohol  and  precipitated  with  ethanol- 
sodium  acetate.  PCR  primers  for  all  834  genes  covering  the  entire  open  reading 
frame  (ORF)  of  each  gene  from  R.  prowazekii  Madrid  E  were  purchased  from 
Research  Genetics  (Huntsville,  Ala.).  PCRs  were  carried  out  in  100-jjlI  reaction 
mixtures  using  a  96-well  PE-9700  thermocycler  (Applied  Biosystems,  Foster  City, 
Calif.).  The  reaction  mixture  contained  15  ng  of  DNA,  0.1  mM  deoxynucleoside 
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triphosphate,  100  ng  of  each  primer,  IX  PCR  buffer,  1.5  mM  MgCl2,  and  2.5  U 
of  Taq  polymerase  (Invitrogen-Life  Technology,  San  Diego,  Calif.).  The  reaction 
conditions  were  as  follows:  initial  denaturation  at  94°C  for  30  s  and  33  cycles  of 
denaturation  (94°C  for  45  s),  annealing  (53°C  for  45  s),  and  extension  (72°C  for 
45  s).  Final  extension  was  performed  at  72°C  for  10  min.  Most  PCR  products 
were  0.5  to  2  kb  in  size.  For  —147  genes  that  were  not  amplified  or  that  yielded 
nonspecific  products  or  faint  bands  under  the  conditions  described  above,  the 
PCR  conditions  were  adjusted  or  new  primers  were  designed. 

(ii)  Preparation  of  glass  slides  for  microarray.  Microscope  slides  were  se¬ 
lected  manually  in  order  to  exclude  those  slides  with  dust  or  scratches.  The  slides 
were  submerged  in  NaOH-ethanol  for  2  h  and  then  coated  with  polylysine  after 
the  residual  NaOH  was  completely  removed  by  vigorously  washing  the  slides  for 
20  min  in  double-distilled  H20.  The  slides  were  immediately  spun  at  550  rpm  for 
3  min  and  then  stored  in  desiccators  at  room  temperature  until  spotting  was 
performed. 

(iii)  Spotting  and  posttreatment  of  slides.  All  PCR  products  were  purified  by 
ethanol-acetate  precipitation  and  robotically  printed  onto  polylysine-coated 
slides  by  the  computerized  OMNI  Grid  Arrayer  (GeneMachines,  San  Carlos, 
Calif.).  The  eight-pin  model  was  used  for  spotting  the  DNAs  of  834  genes  plus 
a  96-gene  duplicate  as  the  internal  control.  The  spotted  DNA  was  cross-linked  to 
polylysine  by  UV  irradiation  and  treated  with  succinic  anhydride  to  minimize 
background  staining. 

Hybridization  of  genomic  DNA.  Genomic  DNAs  from  the  Madrid  E  and 
Breinl  strains  were  used  as  templates  for  direct  incorporation  of  fluorescent 
nucleotide  analogs  (Cy3  and  Cy5  dCTPs;  Amersham,  Piscataway,  N.J.)  by  a 
random-priming  polymerization  reaction  (6).  After  denaturation  of  the  DNA, 
the  two  50-|xl  labeling  reaction  mixtures  containing  2  fig  of  template  DNA  (one 
from  Madrid  E  and  another  from  Breinl);  IX  random-primer  reaction  buffer 
(Invitrogen-Life  Technology);  1.2  mM  (each)  dATP,  dGTP,  and  dTTP;  0.6  mM 
dCTP;  3.75  nmol  of  Cy3-dCTP  or  Cy5-dCTP;  and  40  U  of  Klenow  fragment 
(Invitrogen-Life  Technology)  were  incubated  at  37°C  for  2  h.  By  adding  Cot-1 
DNA,  yeast  tRNA,  and  Tris-EDTA  buffer,  the  labeled  DNA  was  then  concen¬ 
trated,  denatured,  and  applied  to  microarray  slides.  The  hybridization  was  car¬ 
ried  out  at  65°C  overnight,  and  then  the  slides  were  washed  for  5  min  in  2x  SSC 
buffer  (IX  SSC  is  0.15  M  NaCl  plus  0.015  M  sodium  citrate)  with  0.03%  sodium 
dodecyl  sulfate,  followed  by  5  min  in  IX  SSC  and  then  1  to  2  min  in  0.2X  SSC 
buffer.  In  a  subsequent  experiment,  each  DNA  was  reciprocally  labeled  in  order 
to  eliminate  artifacts  due  to  differences  in  labeling  intensity. 

Imaging  and  data  analysis.  The  hybridized  arrays  were  scanned  at  lO-p-m 
resolution  on  a  GenePix  4000  scanner  (Axon  Instruments,  Inc.,  Foster  City, 
Calif.)  at  variable  photomultiplier  tube  voltage  settings  to  obtain  maximal  signal 
intensities  with  <1%  probe  saturation.  The  Cy5-  and  Cy3-labeled  DNA  samples 
were  scanned  at  635  and  532  nm,  respectively.  The  resulting  TIFF  images  were 
analyzed  by  Gene  Pix  Pro  version  3.0.  For  each  spot,  the  background  intensity 
was  subtracted  from  the  spot  signal  intensity  to  produce  a  channel-specific  value. 
The  Cy3/Cy5  ratios  for  all  targets  were  determined  after  background  subtraction. 
Gene  cluster  analysis  was  performed  using  the  National  Cancer  Institute-Center 
for  Information  Technology  microArray  database  system.  Hierarchical  clustering 
was  performed  according  to  the  method  of  Eisen  et  al.  using  an  “uncentered” 
correlation  algorithm  (14).  The  cluster  data  were  filtered  for  signal  intensities 
and  background  of  >  1.5-fold  in  both  channels.  The  raw  data  were  converted  into 
log2  data  for  the  analysis  of  the  genomic  distribution  versus  log2  ratio. 

Southern  blot  analysis.  Southern  blot  hybridization  was  carried  out  by  using 
standard  protocols  (35)  with  the  following  modifications.  Probes  were  prepared 
using  a  PCR  DNA  biotinylation  kit  (Kirkegaard  &  Perry  Laboratories,  Gaith¬ 
ersburg,  Md.)  that  yielded  products  between  0.5  and  1.0  kb  in  size.  MacVector 
version  7.0  software  (Genetics  Computer  Group,  Madison,  Wis.)  was  used  to 
select  restriction  enzymes  for  the  digestion  of  genomic  DNA  to  generate  proper 
length  distributions  and  to  predict  signal  sizes  in  Southern  blot  analysis.  Rick¬ 
ettsial  DNA  was  digested  with  EcoRl/Xbal  and  Hindlll  for  the  detection  of  ftsK, 
ifbE,  IpxA,  and  tpoH.  Hindlll  was  replaced  by  EcoKV /Clal  for  the  detection  of 
virB4.  None  of  these  restriction  enzymes  cut  Madrid  E  DNA  sequences  that  were 
designed  for  the  generation  of  probes.  The  DNA  fragments  were  separated  on 
0.8%  agarose  gels  by  overnight  electrophoresis  and  then  transferred  onto  a 
positive  nylon  membrane  (Kirkegaard  &  Perry  Laboratories)  overnight  by  alka¬ 
line  capillary  transfer.  The  labeled  DNA  probe  was  hybridized  to  the  membrane, 
and  signal  was  detected  with  a  Random  Primer  DNA  biotinylation  kit  (Kirkeg¬ 
aard  &  Perry  Laboratories)  according  to  the  manufacturer’s  recommendations. 

PCR  and  RT-PCR  for  rp084  gene.  PCR  and  RT-PCR  were  performed  with 
primers  specific  for  the  entire  ip084  gene  (forward,  5'  ATG  ACT  ATA  ACA 
GGC  AAA  TT  3';  reverse,  5'  TAT  TTT  TCT  TTT  GCC  TCA  CT  3').  The 
reaction  conditions  were  as  follows:  initial  denaturation  at  94°C  for  30  s;  32  cycles 
of  denaturation  (94°C  for  45  s),  annealing  (52°C  for  45  s),  and  extension  (72°C 
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TABLE  1.  Primers  and  probes"  used  in  real-time  PCR 
Target  Sequence  (5'-3') 


virB4 

Forward . CCT  TTG  CAT  CTT  GAA  GAT  ACT  AGT  GAG  A 

Reverse . TAT  CTT  GCG  CAG  TTA  AAC  TTT  CAC  C 

Probe . AAG  TTA  CAA  GCA  CGC  GAA  GCC  ACT  CTA  AA 


ftsK 

Forward . GCA  GCC  TTA  TTA  TAG  CCG  ATT  CTA 

Reverse . AAA  AGC  AGT  ACA  GAT  AGT  GCG  TAA  TG 

Probe . ACG  TAA  ATC  CTC  CAT  TAG  CTA  TAT  CCA  AAG  GT 

rfbE 

Forward . TCT  TGT  TAT  GAG  GGA  GAA  AAA  ATA  GCT 

Reverse . ATT  CCA  GCA  ATC  AGT  TTT  AG  A  AGT  GA 

Probe . TTT  TAC  CTG  AAC  CAT  TAC  TCC  CAA  TAA 


rpoH 

Forward . TTA  CCT  TCT  GAA  ACC  AAT  TCG  GTA  A 

Reverse . CGA  GTC  ATC  TTA  AAC  TTG  TTG  CAA  A 

Probe . ATT  GCA  AGC  AGC  TAC  AAA  ACT  TAC  GGT  CTT  CCT 


IpxA 

Forward . AAG  CAC  TGC  GCG  TTT  ACT  ACT  TAC 

Reverse . ATT  CGA  TGA  TAG  GCG  GAC  TAT  CA 

Probe . TCC  AAA  CGG  TAT  TAC  ATC  TGC  GCC  AAC  T 


rpl05 

Forward . TCA  GTG  CAT  GAA  CAC  GGT  AAG  TC 

Reverse . GCT  GCG  CTC  CAT  CTT  TAC  G 

Probe . TAC  CAA  TAT  CTG  GCT  TAT  CTC  CCA  AA 


ampG2 

Forward . GCG  ATA  AAC  ATC  CTG  CCA  AGA 

Reverse . GCA  TAA  GTG  AAG  CGG  CTT  GTA  CT 

Probe . AGC  TCC  AGA  ACC  AGA  AAT  AAG  TAT  ACC  AAT  TCT  AAA  ACC  A 


pal 

Forward . AGC  TGC  TGC  CCT  TCT  TTC  G 

Reverse . AGC  TTG  CTG  GTT  ATC  AAA  ACA  TCC 

Probe . AGT  ACC  TCT  TTC  ATC  ACA  ATG  CCC  TTC  AAC  AGT  A 


lies 

Forward . AGG  CGA  CGG  AAC  AGG  AAT  AGT  A 

Reverse . TCA  ACA  GGG  CAA  ACA  AGT  TCA  A 

Probe . TAT  GGC  TCC  AGG  CTT  TGG  TGA  AGA  TGA 


uvrB 

Forward . CAA  CTC  CAC  CTC  AAT  GCT  TGA  TAT 

Reverse . CAG  GTG  AGA  AAC  TTA  CTC  AGC  TTG  ATA 

Probe . CCT  TTC  TCG  TGG  CAT  AAC  AAA  GTG  TGA  ATT  ACC 

16S  RNA 

Forward . GGA  GCA  TGC  GGT  TTA  ATT  CG 

Reverse . AAA  GCA  TCT  CTG  CAA  TCC  GTA  AC 

Probe . ACG  CGA  AAA  ACC  TTA  CCA  ACC  CTT  GAC  AT 


"  Probes  were  sense  or  antisense  sequences  labeled  with  FAM  at  the  5'  end  and  with  TAMRA  at  the  3'  end. 


for  45  s);  and  a  final  extension  step  at  72°C  for  10  min.  Amplification  of  the  16 
S  rRNA  gene  was  done  as  an  internal  control  under  the  same  conditions  as  for 
i-p084  with  the  following  primers:  5'  GTT  CGG  AAT  TAC  TGG  GCG  TA  3' 
(forward)  and  5'  AAT  TAA  ACC  GCA  TGC  TCC  AC  3'  (reverse). 

DNA  sequencing.  PCR  primers  that  flanked  the  sequence  of  Madrid  E  from 
the  end  of  ip082  to  the  beginning  of  tp085  were  designed.  A  PCR  product  from 
Madrid  E  was  calculated  to  be  a  4.3-kb  fragment.  The  PCR  product  from  Breinl  was 

. 2.5  kb  and  was  directly  sequenced  with  the  Big-Dye  terminator  cycle-sequencing 

kit  using  the  automated  sequencing  machine  (Applied  Biosystems).  In  addition,  the 
full-length  genes  virB4,ftsK ,  IpxA,  rfbE,  and  rpoH  of  Breinl  and  Madrid  E  were  PCR 
amplified  and  cloned  into  the  BlueScript  vector  (Stratagene,  La  Jolla,  Calif.)  and 


sequenced  as  described  above.  DNA  sequence  analyses  and  alignments  were  per¬ 
formed  with  MacVector  version  7.0  and  Blast  software. 

RNA  extraction  and  RT.  Total  RNA  was  isolated  from  the  Renografin-purified 
rickettsial  strains  Madrid  E  and  Breinl  with  TRIZOL  (Invitrogen-Life  Technol¬ 
ogy)  according  to  the  manufacturer’s  instructions.  The  RNA  samples  were 
treated  with  DNase  I  (Invitrogen-Life  Technology)  at  37°C  for  1  h  and  then 
extracted  with  phenol-chloroform-isoamyl  alcohol  two  or  three  times.  The  qual¬ 
ity  and  quantity  of  RNA  were  confirmed  by  visualization  of  the  agarose  gel  and 
spectrophotometer  analysis.  Two  micrograms  of  total  RNA  was  used  for  cDNA 
synthesis  at  42°C  for  50  min  with  Superscript  II  (Invitrogen-Life  Technology), 
followed  by  RNase  H  treatment  at  37°C  for  20  min.  A  reaction  with  no  reverse 
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FIG.  1.  DNA  microarray  analysis  of  834  R.  prowazekii  genes  plus  a 
96-gene  repeat  as  internal  control.  (A)  Cohybridization  of  Cy3-labeled 
Madrid  E  DNA  and  Cy5-labeled  Breinl  DNA.  (B)  Reciprocal-labeling 
experiment.  (C  and  D)  Higher  magnifications  of  representative  images 
from  panels  A  and  B.  (E)  Locations  of  all  834  ORFs  within  the  1.1  Mb 
in  Madrid  E  are  mapped  on  the  x  axis.  The  log2  Cy5/Cy3  ratio  for  each 
gene  is  shown  on  they  axis.  Cy5/Cy3  ratios  of  >1.5  or  <0.667,  equiv¬ 
alent  to  ±  0.58  in  the  log2  ratio,  are  indicated  by  the  red  bars. 


ID 

Gene 

Ratio  1 

Ratio2 

RP107 

unknown 

1.817 

1.597 

RP105 

unknown 

2.240 

1.736 

RP303 

rpoH 

1.730 

1.512 

RP103 

virB4 

1.677 

1.675 

RP100 

rpmE 

1 .559 

1.609 

RP007 

IpxA 

1.783 

1.955 

RP868 

glnQ2 

1.561 

1.681 

RP003 

rfbE 

1.587 

1.812 

RP668 

ampG2 

1.760 

1.864 

RP870 

unknown 

1.523 

1.639 

RP825 

mraY2 

1.524 

1.610 

RP725 

unknown 

1.526 

1.585 

RP456 

rpmJ 

1.792 

1.964 

RP617 

ileS 

1.617 

1.702 

RP823 

ftsK 

1.684 

1.616 

RP158 

unknown 

1.690 

1.559 

RP771 

pal 

1.772 

1.618 

RP670 

cspA 

1.891 

1.731 

RP401 

unknown 

1.547 

1.518 

RP451 

sca3 

1.850 

1.568 

RP660 

rpsJ 

1 .905 

1 .552 

RP773 

rpmF 

1.875 

1.594 

RP203 

uvrB 

1.834 

1.557 

RP141 

rpoC 

0.617 

1.702 

RP084 

unknown 

0.158 

0.084 

FIG.  2.  Cluster  image  of  genes  with  fluorescence  ratios  of  ^1.5  or 
<0.667  in  comparisons  of  Madrid  E  and  Breinl  in  two  separate  exper¬ 
iments.  In  lane  1,  Madrid  E  DNA  was  labeled  by  Cy3  (green)  and 
Breinl  DNA  was  labeled  by  Cy5  (red).  Lane  2  shows  a  repeat  experi¬ 
ment  by  reciprocal  labeling.  Ratio  1  and  ratio  2  represent  the  values  of 
25  genes  obtained  from  the  two  individual  experiments. 


transcriptase  was  included  as  a  control.  All  procedures  were  performed  for  both 
strains  under  identical  conditions  to  allow  comparison. 

Real-time  RT-PCR.  The  primer  pairs  and  probes  of  selected  genes  for  real¬ 
time  PCR  were  designed  using  Primer  Express  software  (ABI,  Weiterstadt, 
Germany)  as  shown  in  Table  1.  The  probes  were  labeled  with  6-carboxyfluores- 
cein  and  6-carboxytetramethyl-rhodamine  at  the  5'  and  3'  ends,  respectively. 
PCR  cycling  was  conducted  using  the  SmartCycler  (Cepheid,  Sunnyvale,  Calif.) 
in  reaction  mixtures  with  total  volumes  of  25  p.1  containing  IX  TaqMan  Univer¬ 
sal  PCR  Master  Mix  (Applied  Biosystems),  0.2  p.M  each  primer,  0.1  (xM  probe, 
and  1  fd  of  diluted  cDNA  template.  After  incubation  at  50°C  for  2  min  and 
denaturation  at  95°C  for  10  min,  the  reaction  went  to  45  to  60  cycles  of  dena- 
turation  (94°C  for  15  s)  and  extension  (60°C  for  1  min).  The  16S  rRNA  gene  was 
included  as  an  endogenous  reference,  and  the  comparative  CT  (threshold  cycles) 
method  was  applied  by  using  arithmetic  formulas  (2-AACr)  (ABI  system;  AB, 
Piscataway,  N.J.).  In  this  method,  the  amount  of  target  was  normalized  to  that  of 
the  reference  gene  relative  to  the  calibrator.  The  efficiency  of  the  target  ampli¬ 
fication  and  the  efficiency  of  the  reference  amplification  were  tested  by  assessing 
the  variations  of  A CT  at  different  amounts  of  template  from  the  reference  and 
sample.  No-cDNA  templates  in  PCR  and  RT  products  without  reverse  transcrip¬ 
tase  were  used  as  controls.  Real-time  RT-PCR  assays  were  performed  in  tripli¬ 
cate  for  each  sample,  and  a  mean  value  and  standard  deviation  were  calculated 
for  the  relative  RNA  expression  levels. 
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TABLE  2.  Functional  categories  for  25  genes  identified  by  microarray 


Category" 

Protein  type 

Gene 

Description 

Cell  envelope 

Surface  polysaccharides, 

RP003  rfbE 

O-antigen  ABC  export  system;  ATP- 

lipopolysaccharides,  and 
antigens 

binding  protein 

RP007  IpxA 

UDP-GlcNAc  acyltransferase 

RP771  pal 

Peptidoglycan-associated  lipoprotein 

RP451  sca3 

Cell  surface  antigen 

Murein  sacculus 

RP825  mraY2 

Phospho-IV-acetylmuramoyl-pentapeptide 

transferase 

Transport  and  binding  protein 

RP868  gbiQ2 

Glutamine  ABC  transporter;  ATP- 

binding  protein 

Cellular  processes 

Cell  division 

RP823  ftsK 

Cell  division  protein 

Chaperone  and  stress- 

RP670  cspA 

Cold  shock  protein 

induced 

Translation 

Aminoacyl-tRNA  synthetase 

RP617  lies 

Probably  isoleucyl-tRNA  synthetase 

Ribosomal  proteins;  synthesis 

RP100  ipmE 

Ribosomal  protein  L31 

and  modification 

RP773  rpmF 

Ribosomal  protein  L32 

RP456  ipmJ 

Ribosomal  protein  L36 

RP660  ipsJ 

Ribosomal  protein  S10 

Transcription 

RNA  synthesis  and 

RP141  rpoC 

RNA  polymerase  beta  subunit 

modification 

RP303  rpoH 

RNA  polymerase  sigma-32  factor 

Replication 

RP203  uvrB 

Repair  excision  nuclease  subunit  B 

Other 

RP103  virB4 

Virulence  protein 

RP668  ampG2 

AmpG  protein 

Unknown 

RP084,  RP105,  RP107,  RP158, 

RP401,  RP725,  RP870 

“  Functions  compiled  according  to  Andersson  et  al.  (2). 


RESULTS  AND  DISCUSSION 

Comparison  of  genomic  DNAs  from  Madrid  E  and  Breinl  by 
microarray.  The  cohybridization  of  Cy3-labeled  Madrid  E 
DNA  (green  fluorescence  signal)  and  Cy5-labeled  Breinl  DNA 
(red  fluorescence  signal)  was  applied  on  the  glass  microarray 
slides.  The  resulting  signal  intensities  from  the  two  fluoro- 
phores  were  compared  for  each  gene  in  order  to  detect  the 
genomic  differences,  e.g.,  DNA  duplications  or  deletions,  in 
the  two  strains.  It  should  be  noted  that  no  plasmid  has  been 
found  in  R.  prowazekii.  The  differences  in  signal  should  not  be 
due  to  the  presence,  absence,  or  different  copy  numbers  of  a 
plasmid.  Our  results  demonstrated  that  R.  prowazekii  genes  are 
highly  conserved  between  the  two  strains,  because  most  of 
them  appeared  to  have  yellow  fluorescence  signals  (Fig.  1A  to 
D).  However,  quite  a  few  genes  appeared  to  have  varied  sig¬ 
nals.  Figure  IE  shows  a  distribution  of  log2  ratio  versus 
genomic  location  for  all  genes  studied.  Those  with  a  fluores¬ 
cence  ratio  of  >1.5  or  <0.667  were  selected  for  further  eval¬ 
uation.  Ratios  between  0.67  and  1.5  are  not  reliable  for  pre¬ 
dicting  the  presence  of  differences  in  gene  content  in 
microarray  experiments.  A  total  of  24  genes  showed  >1.5-fold- 
but  <2. 0-fold-higher  hybridization  signals,  and  one  gene 
showed  <0. 667-fold-lower  hybridization  signals,  with  virulent 
Breinl  DNA  (Fig.  2)  than  with  Madrid  E  DNA.  Reciprocal¬ 
labeling  experiments  were  conducted  to  rule  out  the  possibility 
that  the  difference  in  intensity  was  due  to  a  variation  in  Cy3  or 
Cy5  incorporation.  The  ratio  of  fluorescence  for  rpl41  ( rpoC ) 
originally  showed  <0. 667-fold  intensity  but  was  confirmed  to 
be  >1.5  (1.702  and  1.601)  in  two  subsequent  experiments  (Fig. 
2).  Thus,  only  one  gene  showed  <0.667-fold  lower  in  the  ratio 
of  the  Breinl  strain  to  the  Madrid  E  strain.  These  25  genes 


accounted  for  ~3%  (25  of  834)  of  all  genes  detected  on  the 
array.  Eighteen  out  of  25  genes  fell  within  functional  categories 
based  upon  sequence  homology  to  other  bacterial  genes  (2). 
The  remaining  seven  genes  are  unmatched  to  any  homologous 
function  (Table  2)  (2). 

Since  the  fluorescence  ratios  of  24  genes  were  >1.5,  we 
evaluated  the  possibility  of  additional  gene  copies  in  Breinl 
compared  to  Madrid  E  in  light  of  a  report  of  different  copy 
numbers  of  homologous  genes.  There  are  54  paralogous  gene 
families  comprising  147  gene  products  reported  in  the  genome 
of  R.  prowazekii  Madrid  E  (2).  This  observation  may  indicate 
specific  adaptations  of  bacteria  to  their  environments.  The 
presence  of  four  copies  of  SpoT  in  R.  prowazekii  and  Rickettsia 
conorii  might  be  related  to  the  adaptation  of  Rickettsiae  to 
environmental  starvation  in  their  vectors  (25).  Other  possible 
factors  contributing  to  the  increased  fluorescence  signals  may 
be  differences  in  the  number  of  gene  repeat  sequences,  gene 
insertions  in  Breinl,  or  deletion  of  genes  in  Madrid  E. 

DNA  sequencing  and  Southern  blot  analysis.  To  rule  out 
sequence  variation  among  the  25  genes  in  Breinl  and  Madrid 
E,  5  of  them  (ftsk,  rfbE,  IpxA,  rpoH,  and  virB4)  were  cloned  and 
sequenced.  The  results  demonstrated  that  the  sequences  of  all 
five  genes  were  identical.  No  insertions,  duplications,  dele¬ 
tions,  or  mutations  were  seen. 

To  evaluate  the  possibility  of  multiple  copies  or  paralogs  of 
the  disparate  genes  and  to  explain  why  the  hybridization  signal 
ratios  were  between  1.5  and  2.0  with  these  genes,  Southern 
hybridizations  were  performed  for  the  above-mentioned  five 
genes  (Fig.  3).  Total  DNAs  were  isolated  from  the  R. 
prowazekii  Madrid  E  and  Breinl  strains,  respectively,  and  then 
digested  with  restriction  enzymes,  which  were  selected  by 
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FIG.  3.  Comparative  Southern  blot  analyses  of  ftsk,  rfbE,  IpxA, 
rpoH,  and  virB4  for  R.  prowazekii  strains  Breinl  (B)  and  Madrid  E  (E). 
(A)  Hybridization  patterns  with  ftsK,  rfbE,  IpxA,  and  rpoH  probes.  M, 
molecular  size  markers.  (B)  Hybridization  patterns  with  virB4  probe. 
The  arrowheads  indicate  additional  hybridized  bands  in  the  Breinl 
DNA.  Lane  M,  molecular  size  markers. 


MacVector  version  7.0  software  based  on  the  published  Ma¬ 
drid  E  sequence.  The  restriction  enzymes  were  predicted  to  cut 
frequently  within  Madrid  E  DNA  but  not  within  the  regions  for 
probes.  The  Southern  blots  were  hybridized  with  five  DNA 
probes  synthesized  by  PCR.  The  resulting  blots  displayed  one 
strong  band  each  for  ftsK,  ifbE,  IpxA,  and  tpoH  in  all  lanes  of 
Madrid  E  and  Breinl  DNAs  digested  with  EcoRl/Xbal  and 
HiruWW,  which  was  consistent  with  the  predicted  single  copies 
of  these  genes  in  the  Madrid  E  genome  (2).  However,  the  blots 
of  Breinl  DNA  digested  by  EcoRl/Xbal  appeared  to  have  an 
additional  weak  band  for  all  four  genes.  These  additional 
bands  were  not  seen  in  the  Breinl  strain  after  Hindlll  diges¬ 
tions.  The  blots  hybridized  with  the  rpl03  probe  for  the  detec¬ 
tion  of  virB4  displayed  two  major  bands  in  Madrid  E.  One  was 
located  at  the  expected  position  of  3,596  bp  containing  rpl03. 
The  second  band  resulted  from  the  hybridization  of  the  probe 
to  another  virB4  ( rp784 )  and  was  located  at  a  predicted  posi¬ 
tion  of  1,426  bp.  However,  the  blots  from  Breinl  DNA  showed 
a  third  band  at  1.3  kb,  in  addition  to  the  two  similar  bands  seen 
in  Madrid  E  for  virB4.  The  genome  sequence  for  Madrid  E  has 
confirmed  two  paralogs  ( rpl03  and  rp784)  of  virB4,  which 
showed  30.9  and  20.3%  DNA  identities  to  Agrobacterium  tu¬ 
mefactions  virB4  (2).  Further  sequencing  alignments  for  rpl03 
and  rp784  showed  only  49  and  25%  identities  in  DNA  and 
protein,  respectively,  dispersed  over  the  whole  ORF,  but  both 
contain  a  consensus  VirB4  domain.  These  findings  with  virB4 
demonstrated  that  probe  rpl03  could  hybridize  to  rp784  in 
addition  to  itself.  Correspondingly,  the  third  band  in  Breinl 
may  represent  another  hybridization  between  rpl03  and  a 
paralog  of  virB4.  All  these  data  support  the  notion  that  there 
may  be  an  additional  paralog(s)  of  virB4,ftsK,  rfbE,  IpxA,  and 
ipoH  in  the  Breinl  strain  that  is  not  found  in  the  Madrid  E 
strain.  Since  paralog  genes  are  not  identical  in  their  DNA 
sequences,  despite  belonging  to  the  same  functional  family,  it 
is  conceivable  that  the  signals  due  to  paralogous  hybridization 
may  be  not  as  strong  as  those  due  to  homologous  hybridiza¬ 
tion.  We  are  not  sure  why  there  are  no  weak  bands  due  to 
paralogs  in  Breinl  DNA  for  ftsK,  rfbE,  IpxA,  and  rpoH  after 
HindlW  digestion,  as  well  as  for  virB4  after  EcoRV/Clal  diges¬ 
tion.  It  is  possible  that  the  paralogs  have  restriction  cutting 
sites  within  the  probe-complementary  regions,  or  it  may  simply 
be  that  due  to  weak  hybridization  the  signals  were  too  faint  to 
be  seen.  These  data  collectively  suggested  that  an  additional 
paralog(s)  of  the  five  genes  might  exist  in  the  virulent  Breinl 
strain,  providing  an  explanation  for  the  observation  of  gene 
hybridization  ratios  between  1.5  and  2.0  from  microarray  ex¬ 
periments. 

Real-time  RT-PCR  analysis.  To  estimate  the  expression  pro¬ 
files  of  these  genes  in  yolk  sac-grown  R.  prowazekii,  we  con¬ 
ducted  a  comparative  quantification  of  the  gene  expression 
profiles  for  10  of  the  genes  with  detected  differences  in  the 
microarray.  Specific  primers  and  probes  were  designed  for 
each  gene  in  order  to  perform  real-time  PCR.  The  target 
sequences  were  amplified  from  both  Madrid  E  and  Breinl. 
Control  experiments  without  reverse  transcriptase  generated 
weak  signals  with  CT  values  near  60  cycles.  Other  controls 
without  cDNA  templates  were  negative  in  all  experiments.  The 
collected  data  were  analyzed  using  the  standard  comparative 
2~AACt  method.  With  this  method,  the  expression  profiles  of 
the  two  strains  were  compared  using  the  16S  rRNA  gene  as  the 
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TABLE  3.  Detection  of  gene  expression  in  R.  prowazekii  by  real-time  PCR 


Gene 

Madrid  E 

Breinl 

Expression  (Breinl  vs. 
Madrid  E)c 

Cr 

z  ±  SD6 

CT 

x  it  SD 

virB4 

36.65 

36.41  ±  0.21 

34.41 

34.74  ±  0.40 

2.7 

36.26 

34.62 

36.32 

35.19 

ftsk 

36.30 

36.19  ±  0.10 

35.75 

35.47  ±  0.25 

1.4 

36.15 

35.40 

36.11 

35.27 

rfbE 

40.88 

42.00  ±  4.04 

36.98 

40.22  ±  4.08 

2.9 

38.64 

38.87 

46.49 

44.80 

rpoH 

33.29 

33.59  ±  0.39 

32.71 

32.84  ±  0.22 

1.4 

33.45 

32.72 

34.03 

33.10 

IpxA 

35.03 

35.03  ±  0.23 

32.30 

32.69  ±  0.41 

4.2 

35.26 

32.67 

34.81 

33.11 

ipl05 

33.76 

33.65  ±  0.11 

32.96 

32.97  ±  0.08 

1.3 

33.54 

33.06 

33.64 

32.90 

ampG2 

35.11 

35.05  ±  0.17 

31.69 

31.86  ±  0.24 

7.6 

35.18 

31.76 

34.85 

32.14 

pal 

30.04 

30.09  ±  0.35 

27.55 

27.77  ±  0.26 

4.2 

29.77 

27.71 

30.46 

28.05 

iles 

34.60 

34.44  ±  0.30 

32.79 

33.15  ±  0.34 

2.0 

34.09 

33.20 

34.63 

33.47 

uvrB 

39.04 

39.39  ±  0.31 

36.00 

35.83  ±  0.15 

9.9 

39.61 

35.71 

39.52 

35.78 

16S  rRNArf 

18.42 

18.38  ±  0.07 

18.00 

18.12  ±  0.14 

18.30 

18.08 

18.43 

18.28 

“  C,  represents  the  fractional  cycle  number  at  which  a  significant  increase  in  fluorescence  signal  above  a  threshold  signal  (horizontal  zero  line)  can  first  be  detected. 
b x,  average  value. 

c  Using  the  comparative  (AAGr)  method. 
d  Internal  control. 


internal  control.  The  results  showed  that  all  real-time  PCRs 
were  positive  for  the  corresponding  target  sequences  with  CT 
values  from  27  to  42.  The  expression  levels  of  all  10  genes  were 
higher  in  the  virulent  Brienl  strain  than  in  the  avirulent  Madrid 
E  strain,  ranging  from  1.3  to  9.9  times  higher.  The  genes  with 
twofold-higher  or  greater  expression  in  Breinl  were  virB4,  rfbE, 
IpxA,  ampG2,pal,  ileS,  and  uvrB.  These  results  are  summarized 
in  Table  3.  The  real-time  PCR  and  microarray  data  for  these 
10  genes  are  illustrated  in  Fig.  4.  The  expression  ratios  of  IpxA, 
ampG2,  pal,  and  uvrB  were  higher  than  those  of  other  genes 
compared  to  the  DNA  ratios  by  microarray  experiments,  sug¬ 
gesting  different  transcriptional  regulations  in  Breinl  than  in 
Madrid  E. 

Analysis  of  rp084.  One  of  the  intriguing  observations  in  this 
study  was  that  the  rp084  gene  exhibited  only  —1/10  (0.121, 
averaging  from  0.158  and  0.084)  of  the  intensity  in  Breinl  DNA 
as  in  Madrid  E  DNA  (Fig.  1  and  2).  To  verify  the  microarray 
result,  amplification  of  rp084  by  PCR  and  RT-PCR  were  per¬ 
formed  for  the  two  strains.  Neither  the  DNA  nor  the  mRNA  of 
this  gene  was  detected  in  Breinl,  but  both  were  present  in 
Madrid  E  (Fig.  5A  and  B).  PCR  amplification  of  the  region 
flanking  rp083,  rp084,  and  nearby  noncoding  regions  produced 
a  4.3-kb  fragment  in  Madrid  E  (expected  size)  and  a  2.5-kb 


product  in  Breinl.  Sequencing  analysis  of  the  amplicon  from 
Breinl  revealed  that  there  was  a  deletion  of  1,812  bp,  including 
part  of  rp083,  all  of  rp084,  the  noncoding  region  between  rp083 
and  rp084,  and  part  of  the  noncoding  region  after  rp084  (Fig. 
5C).  We  are  not  sure  of  the  significance  of  this  deletion  of 
genomic  DNA  for  the  virulent  Breinl  strain.  The  presence  of 
/ p084  was  detected  by  PCR  amplification  in  another  virulent 
strain  of  R.  prowazekii,  MADII,  but  not  in  the  Ananier  strain 
(data  not  shown).  Both  rp083  and  rp084  belong  to  the  category 
of  unknown  function  in  Madrid  E  (2).  There  were  no  obvious 
similarities  in  either  the  primary  DNA  or  protein  sequence 
after  an  exhaustive  database  search.  Following  a  BLAST  da¬ 
tabase  search  with  the  predicted  mass  (18,232  Da)  and  pi 
(9.52)  for  rp084,  three  mitochondrial  proteins,  CTAG  (cyto¬ 
chrome  c  oxidase  assembly  protein),  CYCJ  (cytochrome  c-type 
biogenesis  protein),  and  OSCP  (ATP  synthase  oligomycin  sen¬ 
sitivity  conferral  protein),  with  the  same  pi  and  very  similar 
molecular  masses  of  —18  kDa  were  identified.  They  show  sim¬ 
ilarities  in  the  predicted  secondary  structures  (data  not  shown) 
and  putative  biochemical  characteristics,  such  as  overall  hydro- 
philicity,  a  high  content  of  basic  amino  acids  (higher  pi  value), 
and  a  predicted  a-helical  secondary  structure,  all  of  which  are 
indicative  of  mitochondrion-related  proteins  (17).  CTAG 
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FIG.  4.  Comparison  of  ratios  of  DNA  signals  by  microarray  and  ratios  of  mRNA  expression  levels  by  real-time  RT-PCR  for  selected  genes  in 
R.  prowazekii  strain  Breinl  versus  strain  Madrid  E. 


found  in  the  rickettsia  Cowdria  ruminantium  was  shown  to 
affect  cytochrome  c  oxidase  synthesis  at  a  terminal  stage  (45). 
CYCJ  might  be  required  for  the  biogenesis  of  c-type  cyto¬ 
chromes  in  Bradyrhizobium  japonicum  and  Rhizobium  etli  (33, 
41).  In  Drosophila  melanogaster  and  Saccharomyces  cerevisiae, 
OSCP  was  implicated  in  proton  conduction  during  ATP  syn¬ 
thesis  (4,  9,  34).  Cytochrome  c  is  believed  to  be  an  important 
factor  in  the  initiation  of  the  apoptotic  cascade.  The  cyto¬ 
chrome  c  released  in  response  to  stress  could  activate  a  caspase 
initiator,  caspase  9,  and  trigger  the  subsequent  activation  of 
downstream  executioner  caspases  in  the  mammalian  host  (22, 
52).  Further  functional  studies  of  rp084  are  necessary  to  dem¬ 
onstrate  the  specific  role  of  this  rickettsial  gene. 

Potential  virulence  genes  identified  in  this  study.  The  virB4, 
ftsK,  IpxA,  and  rfbE  genes  are  known  to  be  virulence  or  viru¬ 
lence-related  genes,  and  their  overrepresentation  in  the  viru¬ 
lent  Breinl  strain  might  indicate  that  they  have  roles  in  the 
virulence  of  R.  prowazekii.  Although  there  are  limited  reports 
of  the  study  of  these  genes  in  Rickettsia,  they  have  been  inves¬ 
tigated  in  Escherichia  coli  and  other  bacteria. 

(i)  virB4.  The  R.  prowazekii  Madrid  E  genome  possesses  two 
virB4  paralogs  ( rpl03  and  rp784)  that  are  homologous  to  the 
virB4  genes  encoding  bacterial  transport  systems,  such  as  T- 
DNA  export  in  A.  tumefaciens,  with  DNA  sequence  identities 
of  30.9  and  20.3%,  respectively  (2).  VirB4  is  believed  to  be  an 
essential  virulence  protein  in  A.  tumefaciens  as  a  cytoplasmic 
membrane  ATPase  (11,  16,  23).  VirB4  also  functions  in  plas¬ 
mid  transport  in  E.  coli,  in  Ptl  toxin  export  in  Bordetella  per¬ 
tussis,  in  the  Cag  system  in  Helicobacter  pylori,  and  in  the 
icm/dot  system  in  Legionella  pneumophila  (reference  12  and 
references  therein).  However,  the  role  of  virB4  in  rickettsiae 
has  not  been  elucidated.  In  addition,  other  virB  family  mem¬ 
bers  have  been  reported  in  R.  prowazekii  Madrid  E,  such  as 
virB8-ll  and  virD4  (2). 


(ii)  ftsK.  FtsK  is  predicted  to  be  involved  in  bacterial  cell 
division.  FtsK  may  possess  two  functions:  the  membrane  do¬ 
main  (FtsKN)  targets  the  protein  to  the  closing  septum  (47), 
and  the  cytoplasmic  portion  (FtsIQ.)  segregates  chromosome 
dimers  during  cell  division  (39,  40,  51).  A  mutation  in  ftsK44 
appeared  to  block  the  late  phase  in  the  cell  division  pathway 
(5).  This  information  suggests  that  the  FtsK  protein  plays  a 
crucial  role  in  cell  division  by  triggering  chromosome  segrega¬ 
tion.  There  has  been  no  report  of  the  study  of  ftsK  in  rickett¬ 
siae. 

(iii)  IpxA  and  rfbE.  IpxA  encodes  UDP-lV-acetylglucosamine 
(UDP-GlcNAc)  acyltransferase,  which  is  required  for  the  first 
step  of  lipid  A  (endotoxin)  biosynthesis  (30).  LpxA  catalyzes 
the  transfer  of  an  R -3' -hydroxy acyl  chain  from  R-3-hydroxyacyl- 
acyl  carrier  protein  to  the  glucosamine  3-OH  of  UDP-GlcNAc 
(1).  The  growth  and  virulence  of  most  gram-negative  bacteria 
require  lipid  A  (27).  It  is  also  a  potent  activator  of  innate 
immunity  in  animal  systems  (30,  43).  The  rfbE  gene  encodes  an 
enzyme  necessary  for  the  biosynthesis  of  the  O-antigen  side 
chain  in  E.  coli  0157:H7  (7).  Rickettsial  lipopolysaccharide 
appears  to  be  an  important  immunogen,  as  well  as  a  virulence 
factor  (36,  38).  Thus,  rfbE  has  a  potential  role  in  the  patho¬ 
genesis  of  R.  prowazekii. 

(iv) .  RpoH.  The  RNA  polymerase  sigma-32  factor  (RpoH) 
has  been  characterized  as  a  primary  component  of  the  heat 
shock  response  in  the  survival  of  E.  coli  (8).  In  response  to  a 
sudden  heat  shock  or  other  stresses,  the  levels  of  RpoH  rise 
transiently  through  enhanced  synthesis  and  protein  stabiliza¬ 
tion.  In  R.  prowazekii,  sigma-73  factor  has  been  overexpressed 
in  E.  coli,  and  it  can  stimulate  transcription  by  promoting  RNA 
polymerase  binding  (3).  Rickettsial  heat  shock  proteins,  in¬ 
cluding  DnaK,  GroEL,  and  GroES,  have  been  detected  by 
Western  blot  analysis  (15),  but  the  regulatory  mechanism  for 
heat  shock  or  stress  is  unknown. 
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(B).  RT-PCR :RP084  and  16  S  rRNA 
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FIG.  5.  PCR  (A)  and  RT-PCR  (B)  amplifications  of  rp084  and  16S 
rRNA  genes  for  Madrid  E  and  Breinl  and  sequencing  analysis.  (C)  Di¬ 
agram  of  genes  around  rp083  and  rp084  in  Breinl  and  Madrid  E  strains. 


In  conclusion,  we  have  successfully  assembled  the  first  ge¬ 
nomewide  DNA  microarray  for  R.  prowazekii  and  used  this 
array  to  compare  genomic  variation  between  the  virulent  Bre¬ 
inl  and  avirulent  Madrid  E  strains.  Differences  in  the  number 
of  paralogs  may  be  present  for  a  few  genes,  such  as  virB4,  ftsK , 
tfbE,  IpxA,  and  rpoH ,  in  Breinl  and  Madrid  E.  The  absence  of 
the  whole  rp084  and  partial  rp083  genes  in  Breinl  is  of  un¬ 
known  relevance  to  rickettsial  virulence.  Our  results  in  this 
report  may  provide  new  insights  into  the  molecular  pathogen¬ 
esis  of  R.  prowazekii ,  such  as  the  differences  in  the  numbers  of 
paralog  genes,  that  may  be  useful  for  the  development  of 
improved  diagnostic  assays  and  vaccines.  Further  investiga¬ 
tions  are  planned  using  this  microarray  to  perform  differential 
expression  profile  analysis  of  rickettsial  mRNAs  from  infected 
mammalian  cells  and  arthropod  vectors. 
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